Spatially resolved single-cell atlas of ascidian endostyle provides insight into the origin of vertebrate pharyngeal organs.

The pharyngeal endoderm, an innovation of deuterostome ancestors, contributes to pharyngeal development by influencing the patterning and differentiation of pharyngeal structures in vertebrates; however, the evolutionary origin of the pharyngeal organs in vertebrates is largely unknown. The endostyle, a distinct pharyngeal organ exclusively present in basal chordates, represents a good model for understanding pharyngeal organ origins. Using Stereo-seq and single-cell RNA sequencing, we constructed aspatially resolved single-cell atlas for the endostyle of the ascidian Styela clava. We determined the cell composition of the hemolymphoid region, which illuminates a mixed ancestral structure for the blood and lymphoid system. In addition, we discovered a cluster of hair cell-like cells in zone 3, which has transcriptomic similarity with the hair cells of the vertebrate acoustico-lateralis system. These findings reshape our understanding of the pharynx of the basal chordate and provide insights into the evolutionary origin of multiplexed pharyngeal organs.

Bivalent activity of super-enhancer RNA LINC02454 controls 3D chromatin structure and regulates glioma sensitivity to temozolomide.

Glioma cell sensitivity to temozolomide (TMZ) is critical for effective treatment and correlates with patient survival, although mechanisms underlying this activity are unclear. Here, we reveal a new mechanism used by glioma cells to modulate TMZ sensitivity via regulation of SORBS2 and DDR1 genes by super-enhancer RNA LINC02454. We report that LINC02454 activity increases glioma cell TMZ sensitivity by maintaining long-range chromatin interactions between SORBS2 and the LINC02454 enhancer. By contrast, LINC02454 activity also decreased glioma cell TMZ sensitivity by promoting DDR1 expression. Our study suggests a bivalent function for super-enhancer RNA LINC02454 in regulating glioma cell sensitivity to TMZ.

Three-Dimensional Gene Regulation Network in Glioblastoma Ferroptosis.

Ferroptosis is an iron-dependent form of cell death, which is reported to be associated with glioma progression and drug sensitivity. Targeting ferroptosis is a potential therapeutic approach for glioma. However, the molecular mechanism of glioma cell ferroptosis is not clear. In this study, we profile the change of 3D chromatin structure in glioblastoma ferroptosis by using HiChIP and study the 3D gene regulation network in glioblastoma ferroptosis. A combination of an analysis of HiChIP and RNA-seq data suggests that change of chromatin loops mediated by 3D chromatin structure regulates gene expressions in glioblastoma ferroptosis. Genes that are regulated by 3D chromatin structures include genes that were reported to function in ferroptosis, like HDM2 and TXNRD1. We propose a new regulatory mechanism governing glioblastoma cell ferroptosis by 3D chromatin structure.

Role of TPEN in Amyloid-ß25-35-Induced Neuronal Damage Correlating with Recovery of Intracellular Zn2+ and Intracellular Ca2+ Overloading.

The overproduction of neurotoxic amyloid-ß (Aß) peptides in the brain is a hallmark of Alzheimer's disease (AD). To determine the role of intracellular zinc ion (iZn2+) dysregulation in mediating Aß-related neurotoxicity, this study aimed to investigate whether N, N, N', N'­tetrakis (2­pyridylmethyl) ethylenediamine (TPEN), a Zn2+­specific chelator, could attenuate Aß25-35­induced neurotoxicity and the underlying mechanism. We used the 3-(4, 5-dimethyl-thiazol-2-yl)-2, 5-diphenyltetrazolium bromide assay to measure the viability of primary hippocampal neurons. We also determined intracellular Zn2+ and Ca2+ concentrations, mitochondrial and lysosomal functions, and intracellular reactive oxygen species (ROS) content in hippocampal neurons using live-cell confocal imaging. We detected L-type voltage-gated calcium channel currents (L-ICa) in hippocampal neurons using the whole­cell patch­clamp technique. Furthermore, we measured the mRNA expression levels of proteins related to the iZn2+ buffer system (ZnT-3, MT-3) and voltage-gated calcium channels (Cav1.2, Cav1.3) in hippocampal neurons using RT-PCR. The results showed that TPEN attenuated Aß25-35­induced neuronal death, relieved the Aß25-35­induced increase in intracellular Zn2+ and Ca2+ concentrations; reversed the Aß25-35­induced increase in ROS content, the Aß25-35­induced increase in the L-ICa peak amplitude at different membrane potentials, the Aß25-35­induced the dysfunction of the mitochondria and lysosomes, and the Aß25-35­induced decrease in ZnT-3 and MT-3 mRNA expressions; and increased the Cav1.2 mRNA expression in the hippocampal neurons. These results suggest that TPEN, the Zn2+-specific chelator, attenuated Aß25-35­induced neuronal damage, correlating with the recovery of intracellular Zn2+ and modulation of abnormal Ca2+-related signaling pathways.

Systematic fine-mapping and functional studies of prostate cancer risk variants.

To date, genome-wide association studies (GWAS) have revealed over 200 genetic risk loci associated with prostate cancer; yet, true disease-causing variants remain elusive. Identification of causal variants and their targets from association signals is complicated by high linkage disequilibrium and limited availability of functional genomics data for specific tissue/cell types. Here, we integrated statistical fine-mapping and functional annotation from prostate-specific epigenomic profiles, 3D genome features, and quantitative trait loci data to distinguish causal variants from associations and identify target genes. Our fine-mapping analysis yielded 3,395 likely causal variants, and multiscale functional annotation linked them to 487 target genes. We prioritized rs10486567 as a genome-wide top-ranked SNP and predicted HOTTIP as its target. Deletion of the rs10486567-associated enhancer in prostate cancer cells decreased their capacity for invasive migration. HOTTIP overexpression in enhancer-KO cell lines rescued defective invasive migration. Furthermore, we found that rs10486567 regulates HOTTIP through allele-specific long-range chromatin interaction.

An autoimmune pleiotropic SNP modulates IRF5 alternative promoter usage through ZBTB3-mediated chromatin looping.

Genetic sharing is extensively observed for autoimmune diseases, but the causal variants and their underlying molecular mechanisms remain largely unknown. Through systematic investigation of autoimmune disease pleiotropic loci, we found most of these shared genetic effects are transmitted from regulatory code. We used an evidence-based strategy to functionally prioritize causal pleiotropic variants and identify their target genes. A top-ranked pleiotropic variant, rs4728142, yielded many lines of evidence as being causal. Mechanistically, the rs4728142-containing region interacts with the IRF5 alternative promoter in an allele-specific manner and orchestrates its upstream enhancer to regulate IRF5 alternative promoter usage through chromatin looping. A putative structural regulator, ZBTB3, mediates the allele-specific loop to promote IRF5-short transcript expression at the rs4728142 risk allele, resulting in IRF5 overactivation and M1 macrophage polarization. Together, our findings establish a causal mechanism between the regulatory variant and fine-scale molecular phenotype underlying the dysfunction of pleiotropic genes in human autoimmunity.

Long range inter-chromosomal interaction of Oct4 distal enhancer loci regulates ESCs pluripotency.

Nuclear architecture underlies the transcriptional programs within the cell to establish cell identity. As previously demonstrated, long-range chromatin interactions of the Oct4 distal enhancer (DE) are correlated with active transcription in naïve state embryonic stem cells. Here, we identify and characterize extreme long-range interactions of the Oct4 DE through a novel CRISPR labeling technique we developed and chromosome conformation capture to identify lethal giant larvae 2 (Llgl2) and growth factor receptor-bound protein 7 (Grb7) as putative functional interacting target genes in different chromosomes. We show that the Oct4 DE directly regulates expression of Llgl2 and Grb7 in addition to Oct4. Expression of Llgl2 and Grb7 closely correlates with the pluripotent state, where knock down of either result in loss of pluripotency, and overexpression enhances somatic cell reprogramming. We demonstrated that biologically important interactions of the Oct4 DE can occur at extreme distances that are necessary for the maintenance of the pluripotent state.

KLF4 inhibits early neural differentiation of ESCs by coordinating specific 3D chromatin structure.

Neural differentiation of embryonic stem cells (ESCs) requires precisely orchestrated gene regulation, a process governed in part by changes in 3D chromatin structure. How these changes regulate gene expression in this context remains unclear. In this study, we observed enrichment of the transcription factor KLF4 at some poised or closed enhancers at TSS-linked regions of genes associated with neural differentiation. Combination analysis of ChIP, HiChIP and RNA-seq data indicated that KLF4 loss in ESCs induced changes in 3D chromatin structure, including increased chromatin interaction loops between neural differentiation-associated genes and active enhancers, leading to upregulated expression of neural differentiation-associated genes and therefore early neural differentiation. This study suggests KLF4 inhibits early neural differentiation by regulation of 3D chromatin structure, which is a new mechanism of early neural differentiation.

Histone chaperone HIRA complex regulates retrotransposons in embryonic stem cells.

BACKGROUND: Histone cell cycle regulator (HIRA) complex is an important histone chaperone that mediates the deposition of the H3.3 histone variant onto chromatin independently from DNA synthesis. However, it is still unknown whether it participates in the expression control of retrotransposons and cell fate determination. METHODS: We screened the role of HIRA complex members in repressing the expression of retrotransposons by shRNA depletion in embryonic stem cells (ESCs) followed by RT-qPCR. RNA-seq was used to study the expression profiles after depletion of individual HIRA member. RT-qPCR and western blot were used to determine overexpression of HIRA complex members. Chromatin immunoprecipitation (ChIP)-qPCR was used to find the binding of H3.3, HIRA members to chromatin. Co-immunoprecipitation was used to identify the interaction between Hira mutant and Ubn2. ChIP-qPCR was used to identify H3.3 deposition change and western blot of chromatin extract was used to validate the epigenetic change. Bioinformatics analysis was applied for the analysis of available ChIP-seq data. RESULTS: We revealed that Hira, Ubn2, and Ubn1 were the main repressors of 2-cell marker retrotransposon MERVL among HIRA complex members. Surprisingly, Ubn2 and Hira targeted different groups of retrotransposons and retrotransposon-derived long noncoding RNAs (lncRNAs), despite that they partially shared target genes. Furthermore, Ubn2 prevented ESCs to gain a 2-cell like state or activate trophectodermal genes upon differentiation. Mechanistically, Ubn2 and Hira suppressed retrotransposons by regulating the deposition of histone H3.3. Decreased H3.3 deposition, that was associated with the loss of Ubn2 or Hira, caused the reduction of H3K9me2 and H3K9me3, which are known repressive marks of retrotransposons. CONCLUSIONS: Overall, our findings shed light on the distinct roles of HIRA complex members in controlling retrotransposons and cell fate conversion in ESCs.

Enhancer architecture-dependent multilayered transcriptional regulation orchestrates RA signaling-induced early lineage differentiation of ESCs.

Signaling pathway-driven target gene transcription is critical for fate determination of embryonic stem cells (ESCs), but enhancer-dependent transcriptional regulation in these processes remains poorly understood. Here, we report enhancer architecture-dependent multilayered transcriptional regulation at the Halr1-Hoxa1 locus that orchestrates retinoic acid (RA) signaling-induced early lineage differentiation of ESCs. We show that both homeobox A1 (Hoxa1) and Hoxa adjacent long non-coding RNA 1 (Halr1) are identified as direct downstream targets of RA signaling and regulated by RARA/RXRA via RA response elements (RAREs). Chromosome conformation capture-based screens indicate that RA signaling promotes enhancer interactions essential for Hoxa1 and Halr1 expression and mesendoderm differentiation of ESCs. Furthermore, the results also show that HOXA1 promotes expression of Halr1 through binding to enhancer; conversely, loss of Halr1 enhances interaction between Hoxa1 chromatin and four distal enhancers but weakens interaction with chromatin inside the HoxA cluster, leading to RA signaling-induced Hoxa1 overactivation and enhanced endoderm differentiation. These findings reveal complex transcriptional regulation involving synergistic regulation by enhancers, transcription factors and lncRNA. This work provides new insight into intrinsic molecular mechanisms underlying ESC fate determination during RA signaling-induced early differentiation.

TPEN attenuates amyloid-ß25-35-induced neuronal damage with changes in the electrophysiological properties of voltage-gated sodium and potassium channels.

To understand the role of intracellular zinc ion (Zn2+) dysregulation in mediating age-related neurodegenerative changes, particularly neurotoxicity resulting from the generation of excessive neurotoxic amyloid-ß (Aß) peptides, this study aimed to investigate whether N, N, N', N'-tetrakis (2-pyridylmethyl) ethylenediamine (TPEN), a Zn2+-specific chelator, could attenuate Aß25-35-induced neurotoxicity and the underlying electrophysiological mechanism. We used the 3-(4, 5-dimethyl-thiazol-2-yl)-2, 5-diphenyltetrazolium bromide assay to measure the viability of hippocampal neurons and performed single-cell confocal imaging to detect the concentration of Zn2+ in these neurons. Furthermore, we used the whole-cell patch-clamp technique to detect the evoked repetitive action potential (APs), the voltage-gated sodium and potassium (K+) channels of primary hippocampal neurons. The analysis showed that TPEN attenuated Aß25-35-induced neuronal death, reversed the Aß25-35-induced increase in intracellular Zn2+ concentration and the frequency of APs, inhibited the increase in the maximum current density of voltage-activated sodium channel currents induced by Aß25-35, relieved the Aß25-35-induced decrease in the peak amplitude of transient outward K+ currents (IA) and outward-delayed rectifier K+ currents (IDR) at different membrane potentials, and suppressed the steady-state activation and inactivation curves of IA shifted toward the hyperpolarization direction caused by Aß25-35. These results suggest that Aß25-35-induced neuronal damage correlated with Zn2+ dysregulation mediated the electrophysiological changes in the voltage-gated sodium and K+ channels. Moreover, Zn2+-specific chelator-TPEN attenuated Aß25-35-induced neuronal damage by recovering the intracellular Zn2+ concentration.

Two Enhancers Regulate HoxB Genes Expression During Retinoic Acid-Induced Early Embryonic Stem Cells Differentiation Through Long-Range Chromatin Interactions.

Homeobox B cluster (HoxB) genes play important roles in retinoic acid (RA)-induced early embryonic stem cells (ESCs) differentiation. Knowledge of regulation network of HoxB is important to further unveil the mechanism of ESCs differentiation. In this study, we identified two enhancers that were activated by RA treatment and 4C data showed long-range interactions between HoxB genes and the two enhancers. CRISPR/Cas9-mediated individual or compound deletion of the two enhancers significantly inhibits HoxB gene expression, and transcriptome analysis revealed that RA-induced early ESCs differentiation was blocked in the enhancer KO cells. We propose new mechanism by which two enhancers regulate HoxB gene expression by different regulation modes during RA-induced early ESCs differentiation through long-range chromatin interactions.

CTCF-binding element regulates ESC differentiation via orchestrating long-range chromatin interaction between enhancers and HoxA.

Proper expression of Homeobox A cluster genes (HoxA) is essential for embryonic stem cell (ESC) differentiation and individual development. However, mechanisms controlling precise spatiotemporal expression of HoxA during early ESC differentiation remain poorly understood. Herein, we identified a functional CTCF-binding element (CBE+47) closest to the 3'-end of HoxA within the same topologically associated domain (TAD) in ESC. CRISPR-Cas9-mediated deletion of CBE+47 significantly upregulated HoxA expression and enhanced early ESC differentiation induced by retinoic acid (RA) relative to wild-type cells. Mechanistic analysis by chromosome conformation capture assay (Capture-C) revealed that CBE+47 deletion decreased interactions between adjacent enhancers, enabling formation of a relatively loose enhancer-enhancer interaction complex (EEIC), which overall increased interactions between that EEIC and central regions of HoxA chromatin. These findings indicate that CBE+47 organizes chromatin interactions between its adjacent enhancers and HoxA. Furthermore, deletion of those adjacent enhancers synergistically inhibited HoxA activation, suggesting that these enhancers serve as an EEIC required for RA-induced HoxA activation. Collectively, these results provide new insight into RA-induced HoxA expression during early ESC differentiation, also highlight precise regulatory roles of the CTCF-binding element in orchestrating high-order chromatin structure.

Branched-chain amino acid aminotransferase-1 regulates self-renewal and pluripotency of mouse embryonic stem cells through Ras signaling.

The developmental plasticity of embryonic stem cells (ESCs) is mainly controlled by well-characterized transcription factors, but additional factors, especially those related to metabolism that modulate this intrinsic program remain elusive. Here, using whole transcriptome analysis, we identified branched-chain amino acid aminotransferase-1(Bcat1) as highly-expressed in mouse ESCs and dramatically down-regulated upon differentiation. Bcat1 deletion impaired pluripotency and self-renewal in mouse ESCs, while Bcat1 overexpression resulted in robust ESC self-renewal and inhibition of differentiation. Whole genome bisulfite sequencing (WGBS) analysis showed that Bcat1 deletion altered whole genome methylation levels and hence gene expression in multiple pathways. Specifically, Bcat1 deletion increased expression of RAS protein activator like 1(Rasal1), leading to inactivation of Ras-Erk/MAPK signaling, while Rasal1 inhibition rescued defects seen in Bcat1 deleted cells. In summary, we demonstrate that Bcat1 is essential for mouse ESC self-renewal and pluripotency and that this effect is mediated by DNA methylation and the Ras signaling pathway.

Epigenetic modulator inhibition overcomes temozolomide chemoresistance and antagonizes tumor recurrence of glioblastoma.

Glioblastoma multiforme (GBM) heterogeneity causes a greater number of deaths than any other brain tumor, despite the availability of alkylating chemotherapy. GBM stem-like cells (GSCs) contribute to GBM complexity and chemoresistance, but it remains challenging to identify and target GSCs or factors that control their activity. Here, we identified a specific GSC subset and show that activity of these cells is positively regulated by stabilization of methyl CpG binding domain 3 (MBD3) protein. MBD3 binds to CK1A and to BTRCP E3 ubiquitin ligase, triggering MBD3 degradation, suggesting that modulating this circuit could antagonize GBM recurrence. Accordingly, xenograft mice treated with the CK1A activator pyrvinium pamoate (Pyr-Pam) showed enhanced MBD3 degradation in cells expressing high levels of O6-methylguanine-DNA methyltransferase (MGMT) and in GSCs, overcoming temozolomide chemoresistance. Pyr-Pam blocked recruitment of MBD3 and the repressive nucleosome remodeling and deacetylase (NuRD) complex to neurogenesis-associated gene loci and increased acetyl-histone H3 activity and GSC differentiation. We conclude that CK1A/BTRCP/MBD3/NuRD signaling modulates GSC activation and malignancy, and that targeting this signaling could suppress GSC proliferation and GBM recurrence.

Folding Keratin Gene Clusters during Skin Regional Specification.

Regional specification is critical for skin development, regeneration, and evolution. The contribution of epigenetics in this process remains unknown. Here, using avian epidermis, we find two major strategies regulate ß-keratin gene clusters. (1) Over the body, macro-regional specificities (scales, feathers, claws, etc.) established by typical enhancers control five subclusters located within the epidermal differentiation complex on chromosome 25; (2) within a feather, micro-regional specificities are orchestrated by temporospatial chromatin looping of the feather ß-keratin gene cluster on chromosome 27. Analyses suggest a three-factor model for regional specification: competence factors (e.g., AP1) make chromatin accessible, regional specifiers (e.g., Zic1) target specific genome regions, and chromatin regulators (e.g., CTCF and SATBs) establish looping configurations. Gene perturbations disrupt morphogenesis and histo-differentiation. This chicken skin paradigm advances our understanding of how regulation of big gene clusters can set up a two-dimensional body surface map.

lncRNA 5430416N02Rik Promotes the Proliferation of Mouse Embryonic Stem Cells by Activating Mid1 Expression through 3D Chromatin Architecture.

Both 3D chromatin architecture and long non-coding RNAs (lncRNAs) play essential roles in pluripotency maintenance. However, whether lncRNAs are involved in organizing 3D chromatin structure remains largely unexplored. We identified 39 lncRNAs bound by Klf4, among which we further revealed the 5430416N02Rik promoter is a chromatin interaction hub. Knockout of the 5430416N02Rik locus reduces the proliferation rate of embryonic stem cells (ESCs). Moreover, deleting both the promoter and the gene body of 5430416N02Rik causes a more severe proliferation defect and has a more profound impact on the transcriptome than deleting the gene body alone. The reduced proliferation of the 5430416N02Rik locus knockout ESCs is mainly due to the downregulation of Mid1, the expression of which requires the inter-chromosomal interaction between Mid1 and 5430416N02Rik loci. In summary, our data demonstrated that the lncRNA 5430416N02Rik gene locus maintains the fast proliferation of ESCs by activating the expression of Mid1 through chromatin interaction.

HOTAIRM1, an enhancer lncRNA, promotes glioma proliferation by regulating long-range chromatin interactions within HOXA cluster genes.

The long noncoding RNA HOTAIRM1 reportedly plays important roles in acute myeloid leukemia, gastric cancer and colorectal cancer. Here, we analyzed potential function of HOTAIRM1 in glioma and asked whether it participates in long-range chromatin interactions. We monitored expression of HOTAIRM1 in glioma tissues and correlated levels with patient survival using the TCGA dataset. HOTAIRM1 was highly expressed in glioma tissue, with high levels associated with shortened patient survival time. We then suppressed HOTAIRM1 activity in the human glioblastoma U251 line using CRISPR-cas9 to knock in a truncating polyA fragment. Reporter analysis of these and control cells confirmed that the HOTAIRM1 locus serves as an active enhancer. We then performed Capture-C analysis to identify target genes of that locus and applied RNA antisense purification to assess chromatin interactions between the HOTAIRM1 locus and HOXA cluster genes. HOTAIRM1 knockdown in glioma cells decreased proliferation and reduced expression of HOXA cluster genes. HOTAIRM1 regulates long-range interactions between the HOTAIRM1 locus and HOXA genes. Our work suggests a new mechanism by which HOTAIRM1 regulates glioma progression by regulating high-order chromatin structure and could suggest novel therapeutic targets to treat an intractable cancer.